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Abstract. We review four mechanisms for forming brown dwarfs; (i) turbulent fragmentation (producing very low-mass 
prestellar cores); (ii) gravitational instabilities in discs; (iii) dynamical ejection of stellar embryos from their placental cores; 
and (iv) photo-erosion of pre-existing cores in HII regions. We argue (a) that these are simply the mechanisms of low-mass 
star formation, and (b) that they are not mutually exclusive. If, as seems possible, all four mechanisms operate in nature, their 
relative importance may eventually be constrained by their ability to reproduce the binary statistics of brown dwarfs, but this 
will require fully 3-D radiative magneto-hydrodynamic simulations. 
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1. Introduction 

The existence of brown dwarfs was first proposed on theo¬ 
retical grounds by Kumar (1963) and by Hayashi & Nakano 
(1963). However, more than three decades then passed be¬ 
fore brown dwarfs were observed unambiguously (Rebolo et 
al., 1995; Nakajima et al., 1995; Oppenheimer et al. 1995). 
Brown dwarfs are now observed routinely, and are estimated 
to be comparable in number with hydrogen-burning stars. It 
is therefore appropriate to ask how brown dwarfs form, and in 
particular to ascertain (a) whether brown dwarfs form in the 
same way as hydrogen-burning stars, and (b) whether there 
is a clear distinction between the mechanisms that produce 
brown dwarfs and those that produce planets. 

In Section 12] we argue that brown dwarfs do form in the 
same way as stars, on the grounds that their statistical prop¬ 
erties (mass function, binary statistics, clustering properties, 
etc.) appear to form a smooth continuum with those of low- 
mass hydrogen-burning stars. We also suggest that under¬ 
standing how brown dwarfs form is the key to answering 
a fundamental anthropic question, namely, what determines 
the lower mass limit for star formation, and thereby the like¬ 
lihood of long-lived stars with habitable zones. In Section 
121 we consider the formation of brown dwarfs by turbulent 
fragmentation, as suggested by Padoan & Nordlund (2002), 
and we address the question of whether an isolated core of 
brown-dwarf mass formed in this way can cool sufficiently 
fast to condense out. In Section |4] we consider the forma- 
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tion of brown dwarfs by gravitational instabilities in discs. 
We stress that only in massive discs, and at large radii, can 
fragments of a disc contract and cool sufficiently fast to con¬ 
dense out; closer in they are likely to bounce and be shredded. 
We also point out that, in a dense proto-cluster, impulsive 
interactions between discs, or between a disc and a naked 
star, should be common, and may be necessary to ensure 
disc fragmentation. In Section |2] we consider the formation 
of brown dwarfs by the ejection mechanism, as suggested by 
Reipurth & Clarke (2001). We point out that the requirements 
for this mechanism to operate are very general, and therefore 
it is likely to occur in nature, although it is probably not the 
only mechanism forming brown dwarfs, given the difficulty 
it has producing close BD-BD binaries. In Sectionj^we con¬ 
sider the formation of brown dwarfs by photo-erosion of pre¬ 
existing cores which are overrun by HII regions, as suggested 
by Hester et al. (1996). We stress that this a very robust mech¬ 
anism, in the sense that it does not require very fine tuning of 
the parameters; but it is also a very inefficient mechanism, in 
the sense that it requires a very massive inital core to form a 
brown-dwarf, and it clearly cannot deliver the brown dwarfs 
in regions like Taurus. In Section0we summarise our review. 

2. Why brown dwarfs appear to form like 
H-burning stars 

We shall assume that brown dwarfs form in the same way 
as hydrogen-burning stars, i.e. on a dynamical timescale, by 
gravitational instability, and with initially uniform elemental 
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composition (reflecting the elemental composition of the in¬ 
terstellar medium out of which they form). Thus, by implica¬ 
tion, we distinguish brown dwarfs from planets, which form 
on a much longer timescale, by the amalgamation of a rocky 
core and - if circumstances allow - the subsequent accretion 
of a gaseous enevelope, resulting in an initially fractionated 
elemental composition with an overall deficit of volatile/light 
elements. If this is the correct way to view the formation 
of brown dwarfs, and we argue below that it is, then brown 
dwarfs should not be distinguished from stars; many stars fail 
to burn helium, and most fail to burn carbon, without forfeit¬ 
ing the right to be called stars. 

The reason for categorising brown dwarfs as stars is that 
the statistical properties of brown dwarfs appear to form a 
continuum with those of low-mass hydrogen-burning stars. 

(i) IMF. The initial mass function is apparently continu¬ 
ous across the hydrogen-burning limit at ^ 0.075 M^. This is 
not surprising, since the processes which determine the mass 
of a star are presumed to occur at relatively low densities 
and temperatures, long before the protostellar material knows 
whether it will reach sufficiently high temperature to burn hy¬ 
drogen before or after reaching sufficiently high density to be 
supported in perpetuity by electron degeneracy pressure. In 
the light of this continuity, it seems perverse to have to speak 
of ‘The IMF for Stars and Brown Dwarfs’ when ‘The Stellar 
IMF’ is already quite long enough. 

(ii) Clustering properties. In clusters, brown dwarfs ap¬ 
pear to be homogeneously mixed with H-burning stars, and 
their kinematics are also essentially indistinguishable. Al¬ 
though they have been searched for - as possible signatures 
of formation by ejection - neither a greater velocity disper¬ 
sion of brown dwarfs in very young clusters, nor a diaspora 
of brown dwarfs around older clusters, has been found. 

(Hi) Binary statistics. Here we have to distinguish at least 
two types of binary system. 

In the first type of binary system, the primary is a Sun-like 
star and the secondary component is a brown dwarf. Amongst 
this type there is a remarkable lack of close systems (the 
Brown Dwarf Desert). However, at larger separations (semi¬ 
major axis a ~ 100 AU), brown-dwarf secondaries are quite 
common. Moreover, the lack of close low-mass secondaries 
is not confined to brown dwarfs. There appears to be a general 
lack of systems with very low mass-ratios, q ~ 0.1. 

In the second type of binary system, the primary is a 
brown dwarf, and therefore the secondary is also a brown 
dwarf (or possibly even a planetary-mass object, if this dis¬ 
tinction must be made, see below). For brown-dwarf pri¬ 
maries, the multiplicity is estimated to be m ~ 30 to 40 %, 
the distribution of semi-major axes peaks at ~ 4 AU 

with a logarithmic dispersion ~ 0.6, and the mean 

mass-ratio is g ~ 0.7. In comparison, G-dwarf primaries are 
estimated to have m ~ 60 %, ~ 30 AU, ^ ~ 1.6, 

and q ~ 0.3. The implication is that, as the primary mass de¬ 
creases, (i) the multiplicity decreases (but only quite slowly), 
(ii) the distribution of semi-major axes shifts to smaller sepa¬ 
rations and becomes narrower (logarithmically), and (iii) the 
distribution of mass ratios shifts towards unity - with these 


trends continuing across the divide between brown dwarfs 
and H-burning stars. 

In fact, the situation is even more complicated than this, 
since there are several systems in which (a) the primary is a 
close binary with Sun-like components (rather than a single 
Sun-like star), and/or (b) the Sun-like primary is orbited at 
large radius (~ 100 AU) by a close BD-BD binary. However, 
the statistics of these systems are limited. 

(iv) Discs, accretion and outflows. Young brown dwarfs 
are observed to have infrared excesses indicative of circum- 
stellar discs, just like young H-burning stars. From their Ha 
emission-line profiles, there is also evidence for ongoing ac¬ 
cretion onto brown dwarfs, and the inferred accretion rates 
form a continuous distribution with those for H-burning stars, 
fitted approxmately by M ~ lO^^Mgyr”^ (M/M^)^. Fi¬ 
nally, the spectra of brown dwarfs also show forbidden emis¬ 
sion lines suggestive of outflows like those from H-burning 
stars, and recently an outflow from a brown dwarf has been 
resolved spatially. Thus, in the details of their circumstellar 
discs, accretion rates and outflows, young brown dwarfs ap¬ 
pear to mimic H-burning stars very closely, and to differ sig¬ 
nificantly only in scale. 

Given this continuity of statistical properties between 
brown dwarfs and H-burning stars, it is probably unhelpful 
to distinguish brown dwarfs from stars, and in the rest of the 
paper we will only use the H-burning limit at ~ 0.075 Mq as 
one of several reference points in the range of stellar masses. 
The D-burning limit at ~ 0.013 Mg therefore falls in the 
same category. We will then define a star as any object form¬ 
ing on a dynamical timescale, by gravitational instability, 
and therefore with uniform interstellar elemental composi¬ 
tion. With this definition, there is the distinct likelihood of 
a small overlap between the mass range of stars and that of 
planets. Given that in the immediate future we are unlikely 
to know too much more than the masses of the lowest-mass 
objects, and certainly not their internal composition, we will 
simply have to accept that there is a grey area in the range 
0.001 to 0.01 Mg which may harbour both stars and planets. 

It follows that understanding how brown dwarfs form is 
important, not just for its own sake, but because it is the same 
as understanding how very low-mass stars form. Thus brown 
dwarf formation is a key part of understanding why most stars 
have masses in the range 0.01 Mg to 100 Mg, and hence why 
there are lots of hospitable stars like the Sun with long-lived 
habitable zones, and enough heavy elements (C, N, O, Si, 
Mg, Al, Fe, etc.) to produce rocky planets and life. The high- 
mass cut-off is probably due to the fact that radiation pressure 
makes it hard to form the highest-mass stars; and the low- 
mass cut-off is probably due to the opacity limit. By studying 
brown dwarf formation we can attempt to confirm and quan¬ 
tify the low-mass cut-off. 

3. Formation by turbulent fragmentation 

The first possibility that we consider is that the processes 
forming prestellar cores create some prestellar cores with 
very low masses. Very low-mass cores must inevitably spawn 
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very low-mass stars, even if they don’t fragment during col¬ 
lapse. This is the formation mechanism that has been ex¬ 
plored by Padoan & Nordlund (2002). By simulating the de¬ 
velopment of interstellar turbulence, they show that a wide 
range of dense structures is formed. If those structures which 
are dense and coherent enough to be gravitationally unstable 
are identified as prestellar cores, they have a mass spectrum 
very similar to the observed stellar IMF. There is support for 
this scheme from the observations of Motte, Andre & Neri 
(1998) who show that the mass function for cores does in¬ 
deed appear to echo the stellar IMF. However, we note (i) 
that the core mass function should relate more closely to the 
system IMF (rather than the stellar IMF), and (ii) that the 
completeness limit of the core mass function does not extend 
to brown-dwarf masses. Moreover, the simulations of Padoan 
& Nordlund do not include gravity, and they use an isother¬ 
mal equation of state. Therefore they do not address the re¬ 
quirement that dynamically contracting cores must be able to 
radiate away at least half the gravitational potential energy 
being released by condensation. 

This requirement is normally referred to as the Opacity 
Limit (but see Masunaga & Inutsuka, 1999), and is presumed 
to determine the minimum mass for star formation. Tradition¬ 
ally, the Opacity Limit has been evaluated on the basis of the 
3-D hierarchical fragmentation picture developed by Hoyle 
(1953). In this picture, a large protocluster cloud becomes 
Jeans unstable and starts to contract. As long as the sound 
speed in the gas remains approximately constant, the increas¬ 
ing density reduces the Jeans mass, and so separate parts of 
the cloud (sub-clouds) become Jeans unstable and can con¬ 
tract independently of one another. This process repeats itself 
recursively, breaking the original cloud up into ever smaller 
and denser sub-sub...sub-clouds, until the gas becomes so 
opaque that it can no longer radiate away the gravitational en¬ 
ergy being released by contraction. At this stage the gas starts 
to heat up, and fragmentation ceases. This yields a minimum 
mass in the range ~ 0.007 M^ to 0.015 M^ (e.g. Low 
& Lynden-Bell, 1976; Rees, 1976; Silk, 1977). 

However, it appears that 3-D hierarchical fragmentation 
does not work. There is no evidence for its occuring in na¬ 
ture, nor does it occur in numerical simulations of star for¬ 
mation. The reason 3-D hierarchical fragmentation does not 
work probably has to do with the fact that the timescale on 
which a fragment condenses out in 3-D is always longer than 
the timescale on which the parent cloud is contracting. There¬ 
fore fragmentation, if it occurs at all, is only temporary, and 
the fragments are then merged by the overall contraction of 
the parent cloud. The only way to avoid this is to start with 
fragments which are widely spaced, but then the rate of accre¬ 
tion onto a fragment is very high, and even if it starts off with 
mass ~ -^jeans’ many times more mas¬ 

sive by the time its contraction becomes non-linear. Thus the 
values for quoted in the previous paragraph are proba¬ 
bly significant underestimates for hierarchical 3-D fragmen¬ 
tation. 

It is therefore appropriate to revisit the question of the 
minimum mass for fragmentation, but now using a model 
which invokes 2-D one-shot fragmentation of a shock- 


compressed layer. We argue that this model is more relevant 
to the contemporary scenario of ‘star formation in a cross¬ 
ing time’ (Elmegreen, 2000), and in particular to the scenario 
simulated by Padoan & Nordlund (2002). In this scenario 
star formation occurs in molecular clouds wherever two - or 
more - turbulent flows of sufficient density collide with suf¬ 
ficient ram pressure to produce a shock-compressed layer out 
of which prestellar cores can condense. The model is ‘2-D’ 
because fragmentation of a shock-compressed layer is in ef¬ 
fect two-dimensional (the motions which initially assemble a 
fragment are largely in the plane of the layer), and it is ‘one- 
shot’ in the sense of not being hierarchical or recursive. 

A shock-compressed layer is contained by the ram pres¬ 
sure of the inflowing gas, and until it fragments it has a rather 
flat density profile. If we consider the simplest case of a head- 
on collision between two streams of equal density, the result¬ 
ing layer fragments at time whilst it is still accumu¬ 

lating, and the fastest growing fragment has mass 


radius (in the plane of the layer) and half-thickness 

-^PRAG (perpendicular to the plane of the layer) given by 

^FRAG 

^ {alGpvfl\ 

(1) 

^FRAG 

~ {a^lG^pvY/\ 

(2) 

^FRAG 

~ {a^IGpvf/^ 

(3) 

'^FRAG 

~ {a^/Gpv'^f/'^ . 

(4) 

Here a 

is the net velocity dispersion in the layer, p is 

the 


pre-shock density in the colliding flows, and v is the relative 
speed with which the flows collide. We note (a) that the frag¬ 
ments are initially flattened objects (’"prag/'^^prag ^ 

1); (b) that is not simply the standard 3-D Jeans mass 

evaluated at the post-shock density and velocity dispersion - 
it is larger by a factor and (c) that our analysis ig¬ 

nores magnetic fields and the possibility that the post-shock 
gas is turbulent. If present, both magnetic fields and turbu¬ 
lence will act to increase the minimum fragment mass. 

2-D one-shot fragmentation has the advantage that the 
fastest-condensing fragment has finite size, i.e. fragments 
with initial radius ~ fpRAc condense out faster than either 
larger or smaller fragments. Moreover, we can analyse the 
growth of a fragment in a shock-compressed layer, taking 
account of the continuing inflow of matter into the frag¬ 
ment. Hence we can identify the smallest fragment which 
can cool radiatively fast enough to dispose of both the PdV 
work being done by compression, and the energy being dis¬ 
sipated at the accretion shock where matter continues to flow 
into the fragment; these two sources of heat turn out to be 
comparable. We find (Boyd &. Whitworth, 2005) that for 
shocked gas with temperature T ~ 10 K and no turbulence 
(i.e. velocity dispersion, ct, equal to the isothermal sound 
speed, 0.2kms“^), the smallest fragment which can con¬ 
dense out is less than 0.003M^, and fragments with mass 
below 0.005 M^ condense out for a wide range of pre-shock 
density p and shock-speed v (as illustrated on Fig. 1). We em¬ 
phasise that this analysis is more robust than the standard one 
based on 3-D hierarchical fragmentation, on two counts, (i) 
The fragments have condensation timescales shorter than all 
competing length scales (a well known property of layer frag¬ 
mentation, e.g. Larson, 1985), so they do not tend to merge 
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Fig. 1. A log/log plot of the (p, v) plane. The dots mark combinations of pre-shock density, p, and collision speed, v, for 
which the fastest growing fragment has a mass less than 0.005 Mq; we assume that the effective post-shock sound speed 
is cr = 0.2kms“^, corresponding to molecular gas at 10 K. The irregularities in the boundaries of this region reflect the 
tendency of marginally unstable low-mass fragments to undergo pulsations before they collapse. The solid line is the locus 
below which p must fall if our treatment of the radiation from the accretion shock is to be valid; see Boyd & Whitworth 
(2005) for details. 


with their neighbours, (ii) Ongoing accretion is taken into 
account. Indeed, the smallest fragment of all starts off with 
mass 0.0011 M^, and grows to 0.0027 before its contrac¬ 
tion becomes non-linear. We conclude that stars with masses 
down to 0.003 Mg can condense out of shock-compressed 
layers. If the temperature of the post-shock gas can be re¬ 
duced further still, to below ~ 6 K, then it is even possible to 
form ‘stars’ with masses below 0.001 Mg. 

4. Formation by disc fragmentation 

Another possibility is that an initially massive prestellar core 
(i.e. significantly more massive than a brown dwarf) spawns 
brown dwarfs by fragmentation. The fragmentation of col¬ 
lapsing cores is a large and complicated topic. However, one 
of the main fragmentation mechanisms which operates in nu¬ 
merical simulations is that a relatively massive primary pro¬ 
tostar forms, surrounded by a massive disc-like structure (al¬ 
beit not necessarily a relaxed rotationally supported disc). 


and then lower-mass secondary protostars - including proto¬ 
brown-dwarfs - condense out of the disc (e.g. Bate, Bonnell 
& Bromm, 2002a,b, 2003; Hennebelle et ak, 2004; Goodwin 
et al. 2004a,b,c). This is the mechanism of core fragmentation 
on which we shall concentrate here. 

If we consider a relaxed massive disc in isolation, there 
is some doubt as to whether it will fragment gravitationally, 
spawning low-mass companions to the central primary proto¬ 
star, or whether spiral modes will act to quickly redistribute 
angular momentum, thereby stabilising - and ultimately dis¬ 
sipating - the disc before it can fragment. However, if a mas¬ 
sive protostellar disc interacts impulsively with another disc, 
or with a naked star, or if the disc simply never has time to 
relax towards an equilibrium state, then it can be launched 
directly into the non-linear regime of gravitational instabil¬ 
ity, and fragmentation is then much more likely. In the dense 
proto-cluster environment where most protostars are born, 
such impulsive interactions must be quite frequent. Therefore 
Boffin et al. (1998) and Watkins et al. (1998a,b) have simu- 
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lated parabolic interactions between two protostellar discs, 
and between a single protostellar disc and a naked proto¬ 
star. All possible mutual orientations of spin and orbit are 
sampled. The critical parameter turns out to be the effective 
shear viscosity in the disc. If the Shakura-Sunyaev param¬ 
eter is low, ctgg ~ 10“^, most of the secondary protostars 
have masses in the range 0.001 to 0.01 M^. Conversely, 
if ags is larger, agg ~ 10“^, most of the secondary protostars 
have masses in the range 0.01 Mg to 0.1 Mg. The formation 
of low-mass companions is most efficient for interactions in 
which the orbital and spin angular momenta are all parallel; 
on average 2.4 low-mass companions are formed per inter¬ 
action in this case. If the orbital and spin angular momenta 
are randomly orientated with respect to each other, then on 
average 1.2 companions are formed per interaction. 

In the above simulations the gas is assumed to behave 
isothermally, which is probably a reasonable assumption, 
since the discs are large (initial radius 1000 AU) and most 
of the secondary protostars form at large distance (perias- 
tra ~ 100 AU). However, disc fragmentation is probably 
not possible at smaller radii because the ambient tempera¬ 
ture close to the central primary protostar is higher and the 
surface-density of the disc is also higher. Consequently the 
optical depths through proto-fragments are higher and they 
are unable to cool radiatively sufficiently fast to condense out 
(Rafikov, 2005); instead they contract adiabatically, bounce, 
and are shredded by tidal forces. Thus gravitational fragmen¬ 
tation is probably limited to the outer parts of such discs. Rice 
et al. (2003) present SPH simulations of discs fragmentating 
gravitationally at small radii (~ 10 AU), but they use a phe¬ 
nomenological cooling law of the form du/dt = —u/fcooL’ 
and the values of fcooL which they invoke are unrealistically 
short; also their cooling law seems to admit indefinite cool¬ 
ing and their discs appear to fragment only after the cooling 
becomes catastrophic. Boss (2001, 2003) also presents sim¬ 
ulations of discs fragmentating gravitationally at small radii, 
performed using a finite difference code with radiation trans¬ 
port. However, the reality of the fragments he finds is ques¬ 
tionable on two counts. First, in evaluating the boundedness 
of the fragments he appears to neglect their internal kinetic 
energy; in a fragment which is bouncing, or contracting but 
destined to bounce, this can be a dominant term in the Virial 
Theorem. Second, he argues that his fragments are cooling 
by convection, but convection cannot contribute to the cool¬ 
ing of a fragment which is condensing out on a dynamical 
timescale. The velocity fields which Boss attributes to con¬ 
vective motions may actually be due to bouncing - in which 
case they will lead to dissolution of his fragments by shred¬ 
ding. 

This may help to explain the Brown Dwarf Desert. Brown 
dwarf companions to solar-type primaries can form by disc 
fragmentation, but only at large radii. To end up in closer or¬ 
bits, they must either accrete material with low specific an¬ 
gular momentum, which will tend to increase their mass; or 
they must interact dynamically with a third star, but this tends 
to place the more massive star in the close orbit, and to eject 
the less massive star (i.e. the brown dwarf). 


We note that the two mechanisms discussed thus far are 
not mutually exclusive. If the initial prestellar core is already 
of very low mass, then it will inevitably produce a low-mass 
protostar, irrespective of whether it fragments or not. If it has 
higher mass, it can only produce a very low-mass protostar by 
fragmenting, and one possible mode of fragmention involves 
the formation of a disc. 

5. Formation by ejection 

The collapse of a prestellar core is unlikely to lead to a sin¬ 
gle star. Even quite modest levels of turbulence (e.g. Good¬ 
win, Whitworth & Ward-Thompson, 2004a) and/or global 
rotation (Cha & Whitworth, 2003; Hennebelle et al., 2004) 
are sufficient to ensure fragmentation. Hence prestellar cores 
usually spawn small-W clusters of protostars (N ~ 2 to 6; 
e.g. Hubber & Whitworth, 2005), which then grow by com¬ 
petitive accretion and interact dynamically (Whitworth et al., 
1995; Bonnell et al., 2001). Protostars which get ejected from 
the core before they have time to grow to 0.075 Mg end up 
as brown dwarfs (Reipurth & Clarke, 2001). It seems in¬ 
escapable that this mechanism occurs in nature, since all that 
is required is the formation and coexistence of more than 
two protostars in a core, with one of them being less mas¬ 
sive than 0.075 Mg; W-dody dynamics will then almost in¬ 
evitably eject one of the protostars, and usually the least mas¬ 
sive one. 

Several numerical simulations have been performed, us¬ 
ing SPH with sink particles, to demonstrate the viability of 
this mechanism, both in cores with high levels of turbu¬ 
lence (Bate, Bonnell & Bromm, 2002a,b, 2003; Delgado Do¬ 
nate, Clarke & Bate, 2003, 2004), and in cores with low lev¬ 
els of turbulence (Goodwin, Whitworth & Ward-Thompson, 
2004a,c). Many of the brown dwarfs formed in these sim¬ 
ulations retain low-mass discs (M^g^, ~ 0.010 Mg and 

i?Disc ~ 40 AU) even after ejection, from which they con¬ 
tinue to accrete. They also have a radial velocity distribution 
which is scarcely distinguishable from that of the hydrogen¬ 
burning stars. This is firstly because part of the overall veloc¬ 
ity dispersion is due to the motions of the different cores rel¬ 
ative to one another, and this part is inherited by all stars; and 
secondly because the brown dwarfs are ejected with rather 
modest velocities (~ lkms“^), and the higher-mass stars 
involved in the ejection also have increased velocity disper¬ 
sion due to their recoil and their now being in a harder bi¬ 
nary system. The main concern with these simulations is that, 
by invoking sink particles, protostellar embryos are instanta¬ 
neously converted into point masses. This predisposes them 
to dynamical ejection, and prohibits them from merging or 
fragmenting further. Therefore the efficiency of the mecha¬ 
nism may have been overestimated. 

Additional support for the mechanism comes from Good¬ 
win et al. (2004b), who present an ensemble of simulations of 
the collapse and fragmentation of cores having a mass spec¬ 
trum, density profiles, and low levels of turbulence, matched 
to those observed in Taurus. These simulations reproduce 
rather well the unusual stellar IMF observed in Taurus, in- 
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eluding the relative paucity of brown dwarfs. As far as we are 
aware, these are the first simulations to demonstrate a direct 
causal link between the core mass spectrum and the stellar 
IMF. 

Again we note that this mechanism does not exclude the 
previous two; indeed it requires formation by fragmentation 
of a collpasing core as a precursor to produce the low-mass 
protostellar embryos which then get ejected. However, ejec¬ 
tion is unlikely to be the only mechanism forming brown 
dwarfs, since it seems very unlikely to produce the rather 
large numbers of close BD-BD binaries observed. 

6. Formation by photo-erosion 


7. Conclusions 

We have discussed four possible mechanisms for forming 
brown dwarfs: turbulent fragmentation, disc fragmentation, 
dynamical ejection and photo-erosion. None of these is mu¬ 
tually exclusive, and in fact the first three may occur con¬ 
secutively. We emphasise that none of these mechanisms has 
been modelled properly with a fully radiative 3-D magneto- 
hydrodynamical code. Therefore, neither the thermal effects 
which presumably determine the minimum mass for star for¬ 
mation, nor the angular momentum transport processes which 
presumably determine the binary statistics, nor the A^-body 
dynamics which presumably determine the clustering prop¬ 
erties of brown dwarfs, has yet been properly evaluated. 


A fourth - and somewhat separate - mechanism for forming 
brown dwarfs is to start with a pre-existing core of standard 
mass (i.e. ~ )and have it overrun by an HIl region (Hes¬ 
ter et ah, 1996). As a result, an ionisation front (IF) starts to 
eat into the core, ‘photo-eroding’ it. The IF is preceded by a 
compression wave (CW), and when the CW reaches the cen¬ 
tre, a protostar is created, which then grows by accretion. At 
the same time, an expansion wave (EW) is reflected and prop¬ 
agates outwards, setting up the inflow which feeds accretion 
onto the central protostar. The outward propagating EW soon 
meets the inward propagating IE, and shortly thereafter the IE 
finds itself ionising gas which is so tightly bound to the pro¬ 
tostar that it cannot be unbound by the act of ionisation. All 
the material interior to the IE at this juncture ends up in the 
protostar. On the basis of a simple semi-analytic treatment, 
Whitworth & Zinnecker (2004) show that the final mass is 
given by 


M - O.OIM^ ( - 

® V0.3kms-W 


A/;., 


-1/3 


1050 g-1 


V10^cm“3/ ’ 


(5) 


where is the isothermal sound speed in the neutral gas of 
the core, A/j^yc which the star(s) exciting the HII 

region emit hydrogen-ionising photons, and is the density 
in the ambient HII region. 

This mechanism is rather robust, in the sense that it pro¬ 
duces very low-mass stars for a wide range of initial condi¬ 
tions, and these conditions are likely to be realized in nature. 
Indeed, the evaporating gaseous globules (EGGs) identified 
in M16 by Hester et al. (1996) - and subsequently in other 
HII regions - would appear to be pre-existing cores being 
photo-eroded in the manner we have described. However, the 
mechanism is also very inefficient, in the sense that it usu¬ 
ally takes a rather massive pre-existing prestellar core to form 
a single very low-mass star. Moreover, the mechanism can 
only work in the immediate vicinity of an OB star, so it can¬ 
not explain the formation of all brown dwarfs, and another 
mechanism is required to explain those seen in star formation 
regions like Taurus. Brown dwarfs formed in this way should 
include close BD-BD binaries, but they are unlikely to retain 
large accretion discs. 
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